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Abstract 

We consider two linearly coupled masses, where one mass can have in- 
elastic impacts with a fixed, rigid stop. This leads to the study of a two 
degree of freedom, piecewise linear, frictionless, unforced, constrained me- 
chanical system. The system is governed by three types of dynamics: cou- 
pled harmonic oscillation, simple harmonic motion and discrete rebounds. 
Energy is dissipated discontinuously in discrete amounts, through impacts 
with the stop. We prove the existence of a nonzero measure set of orbits 
that lead to infinite impacts with the stop in a finite time. We show how 
to modify the mathematical model so that forward existence and unique- 
ness of solutions for all time is guaranteed. Existence of hybrid periodic 
orbits is shown. A geometrical interpretation of the dynamics based on 
action coordinates is used to visualize numerical simulation results for the 
asymptotic dynamics. 



1 Introduction 



Mechanical impact oscillators are ideal case studies for the investigation of the 
dynamical behaviors possible in systems governed by a mixture of continuous 
and discrete dynamics. Guckenheimer and Holmes [6] (pp. 104-116) discuss the 
dynamics of a ball bouncing on a sinusoidally forced vibrating table. In [11] 
Salapaka et al. study a variation of this system, where an unforced harmonic 
oscillator interacts with a sinusoidally vibrating table. A one degree of freedom 
periodically forced oscillator with a stop restricting its motion is studied by Shaw 
and Holmes [13]. Hindmarsh and Jefferies [7] analyze the impact oscillator for the 
case in which the position of the wall is offset from the rest point of the oscillator. 
New types of bifurcations that arise in connection with grazing impacts have been 
extensively studied; see Whiston [15], Nordmark [10], Chin et al. [2] and Budd 
and Dux [1]. For a study of a two degrees of freedom impact oscillator, see the 
work by J. Shaw and S. Shaw [12]. In [4], Engleder et al. investigate, numerically 
and experimentally, a non-smooth system consisting of two adjacent beams that 
can impact each other, vibrate independently or vibrate as a single unit. The 
paper by Stewart [14] provides a good survey of fundamental problems arising 
in the modelling of hybrid systems. We refer the reader interested in a more 
in depth mathematical discussion on the modelling of discontinuous differential 
equations to Deimling [3] , Fihppov [5] or Monteiro [9] . 

In this paper we analyze the dynamics of the two degrees of freedom coupled 
oscillator constrained by a stop shown in Figure 1. There is no external forcing. 
The springs are linear and frictionless. Mass one and mass two interact only via 
the spring ki2, as in a classical coupled oscillator. The sole impacts in the system 
are therefore those of mass one with the stop. CoUisions with the stop are treated 
using a coefficient of restitution model, i.e., the velocity of the mass just after an 
impact is assumed to be always a constant fraction e of its velocity just before 
the impact. Collisions are assumed to be inelastic, with < e < 1. The stop is 
offset a positive distance away from the equilibrium position of mass one. 
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Figure 1: The two degree of freedom impact oscillator. 
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This is an instance of a hybrid dynamical system, also known as a non-smooth 
dynamical system. There exist three distinct types of motion. If mass one is not 
touching the stop then we have the usual two coupled harmonic oscillators. If 
mass one hits the stop, then there is an instantaneous rebound, governed by the 
coefficient of restitution model. Thirdly, mass one could be constrained to remain 
against the stop by the force acting on it, leaving mass two moving as a single 
harmonic oscillator. 

An interesting feature of the dynamics is the presence of orbits that lead to 
infinite impacts of mass one with the stop in a finite time. We prove the existence 
of such orbits and we show that they can be continued from the limit point 
they were approaching, in order to guarantee existence of solutions for all time. 
Another interesting feature of this system is the way energy is dissipated. Because 
the springs are frictionless and the coefficient of restitution is less than unity, 
energy is dissipated only at impacts of mass one with the stop. Furthermore, at 
such an impact, energy is instantaneously reduced by a discrete amount, since 
the speed of the impacting mass is reduced by a discrete amount. 

The paper is divided into two parts. Modelling and Dynamics of the System. 
The Modelling part is concerned with rigorously establishing forward existence 
and uniqueness of solutions for all time. We start by giving an exact description 
of the system and its governing dynamics. In section 2.2 we prove the existence 
of a nonzero measure set of orbits that lead to infinite rebounds in a finite time. 
Then, in section 2.3, we show how to modify our mathematical model to ensure 
existence and uniqueness of solutions for all time. These two sections are the 
more mathematical ones of the paper. Readers whose interest lies more in the 
general dynamical behavior of our system may consider skipping them. Con- 
versely, readers whose interest is the problem posed by infinite rebounds and how 
to possibly deal with that problem, may consider skipping the last part. Dynam- 
ics of the System, as this part concentrates on the dynamics of our particular 
system. In the Dynamics of the System part, we present some of the more in- 
teresting dynamical behaviors of the system. Section 3.1 is devoted to hybrid 
periodic orbits, a surprising form of periodic behavior where the system alter- 
nates between coupled harmonic motion and simple harmonic motion with mass 
one pressed against the stop. Section 3.2 gives a geometrical interpretation of the 
system dynamics based on action coordinates. Using it, we also visualize some 
numerical simulation results. 

2 Modelling 

2.1 The basic model 

Figure 1 shows the system under study and the coordinate system used. The 
origin of the Xi and X2 coordinates in Figure 1 correspond, respectively, to the 
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equilibrium positions of masses one and two if the stop were not present. Con- 
stants k and ki2 denote respectively the outer and middle spring constants as 
shown in the figure. The forces on mass one and mass two are denoted by Fi and 
F2. For simplicity, we let both masses be unity throughout the paper. The posi- 
tion of the stop, denoted by s, is measured with respect to the reference defined 
by Xi — 0. At the equilibrium position mass one is not leaning against the stop, 
that is, s is assumed to be strictly positive. Also, in all cases it is assumed that 
mass one is located to the left of the stop. 

Depending on the value of the state variables, (xi, fi, 0:2, V2), where Vi = Xi 
and V2 = X2, a different set of equations governs the evolution of the motion. 
Namely, at any time instant one of the following three types of dynamics holds: 

1. Coupled Oscillator Dynamics (COD): The coupled harmonic oscillator 
dynamics governs the system when mass one is not against the stop. 

2. Rebound Map (R): When mass one hits the stop with a nonzero velocity, 
an instantaneous rebound occurs. The rebound is modeled using a coefficient of 
restitution model. The coefficient of restitution is assumed to be strictly between 
zero and one. 

3. Contact Dynamics ( CD): If at some instant mass one reaches the stop with 
zero velocity, then the force on it may constrain it to remain stuck against the 
stop. During such a period, the second mass undergoes simple harmonic motion. 

Formally, we divide the phase space into three regions and assign to each 
region one of the aforementioned dynamics. This is laid out in chart 2.1. 

To simplify the way we denote regions in the phase space, we will use an infor- 
mal notation for sets. For example {s = 10, X2 < c} stands for {x — [xi, Vi, X2, V2) 
& S : Xi = 10 A X2 < c} , where S is the phase space. 



Chart 2.1. 

Phase Space: S = Vi, 0:2, i'2) £ | < s} 

1. Coupled Oscillator Dynamics (COD) 

Region: {xi<s,Vi,X2,V2}U{s,Vi<0,X2,V2}U{s,Vi>0,X2,V2}U{s,0,X2 < 

X2,V2}, 

where X2 — ( ^^f" ) s 

Entrance Surface: {s, Vi < 0, X2, V2} U {s, 0, X2, V2 < 0} 
Exit Surface: {s, Vi > 0, X2, V2} U {s, 0, £2, V2 > 0} 



Dynamics: < 



Xl = Vi 

Vl = -{k + ki2)Xi + ki2X2 
X2 = V2 
. V2^ -{k + ki2)x2 + ki2Xi 
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2. Contact Dynamics (CD) 

Region: {s, 0, X2 > X2, V2} U {s, 0, X2, > 0} U {s, 0, X2, V2 < 0} 
Entrance Surface: {s,0,X2,V2 > 0} 
Exit Surface: {s, 0, X2, V2 < 0} 
' xi = 

X2 = ■J^2 

V2 = -(/S + fcl2)(a;2 - X2eq) where, X2eg = (felfe) « 



Dynamics: < 



3. Rebound Map (R) 

Region: {s, vi > 0, X2, f2} U {s, vi < 0, X2, V2} 
Entrance Surface: {s,Vi > 0,a:2,f2} 
Exit Surface: {5,^1 < 0,^2,^2} 




Xi 



Dynamics: < _ 
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In the above chart, the Entrance and Exit surfaces of a region are, respectively, 
the sets through which orbits enter and leave that region. For the rebound map, 
one defines the entrance surface as the domain of the map and the exit surface 
as its range. 

Let us see how the division of the phase space relates to the actual physical 
dynamics of our impact oscillator. Clearly the region {xi < s} is governed by 
COD. More interesting are the dynamics and division of the phase space at the 
hypersurface {xi = s} which physically corresponds to mass one being in contact 
with the stop. This 3-D hypersurface is shown in Figure 2. 

The {xi = s,Vi> 0} region is the rebound entrance set. After a rebound, Vi 
becomes negative and mass one leaves the stop, hence points in {xi = s,vi< 0} 
fall again under COD. Points on the plane {xi — s, vi — 0} are reached when mass 
one decelerates towards the stop and reaches it with zero speed. In this case, if 
the force, Fi. on mass one is positive then mass one stays pressed against the 
stop, otherwise it does not. The critical position of mass two that yields Fi = 
when mass one is at the stop is X2 = X2. The point X2eq is the equilibrium position 
of mass two when mass one is against the stop. Hence, it is the center of the 
simple harmonic motion of mass two under CD. The continuous lines on Figure 2 
represent these CD orbit segments. Observe that the existence of the CD orbits 
moving on the {xi = s} surface precludes the use of this surface as a Poincare 
cross section where one could define an implicit composite map to describe the 
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Figure 2: The x\ — s hypersurface. 



global dynamics. 

Now, one would like to formally establish forward existence and uniqueness 
of solutions. It turns out that each of our three types of dynamics guaranteeing 
existence and uniqueness within their assigned regions is not enough to establish 
global existence and uniqueness of solutions for all time. In non-smooth systems, 
forward existence can be broken by orbits that lead to an infinite number of 
switches between different types of dynamics in a finite time, for these orbits are 
clearly undefined for times posterior to the occurence of those infinite switches. 
Surprisingly, such orbits do occur in our system and in fact their set has nonzero 
measure. These are orbits that undergo infinite rebounds in a finite time, in other 
words, orbits that switch an infinite number of times between COD dynamics and 
R dynamics in a finite time. In the next section we will prove the existence of 
such orbits and that their set has nonzero measure. Then, in section 2.3, we solve 
the existence problem through continuation of these orbits from the limit point 
they were approaching. We also rule out any other type of infinite switching in 
finite time, hence guaranteeing forward existence and uniqueness of orbits in the 
system. 

2.2 Proof of the existence of a nonzero measure set of 
orbits that lead to infinite rebounds in a finite time 

This section is devoted to proving the theorem below. 

Theorem 1 There exists a nonzero measure set of orbits that lead to infinite 
rebounds in a finite time. 

For reference, we group below all definitions and basic assumptions used in 
this section. 
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Fi{t) = net force on mass one. 
F2{t) = net force on mass two. 

X2 = (^fc^) the position of mass two that results in Fi = when Xi — s. 

X2eq = (fc^^) ■S) the position of mass two that results in F2 = when xi = s. 
Voutn = velocity of mass one just after the nth impact. 
Vin^ = velocity of mass one just before the nth impact. 
Atn = time interval between impacts n and n+1. 
Assume s > and < £ < 1. 

As the proof of Theorem 1 will be constructive, we select an initial condition 
(a;i(0),i;i(0),X2(0),t;2(0)) = (xio, t^io, 2^20, ^^20) such that 

Xio = s 

< (1) 

X2Q > X2 
. ^^20 < 

Also, select a point X2 such that £2 < ^2 < 2:20 • 
Note that 

< X2eq <X2<X2< X20 (2) 

Before turning to the actual proof of Theorem 1, we must prove a series of 
lemmas. We start by proving that lemmas 2 through 8 are true while X2{t) > X2 
holds. 

Lemma 2 < < Fi{t) when the system is under COD, where F^"'" := 

— {k + ki2)s + ki2X2. 

Proof: Under COD, the configuration of mass one and mass two that min- 
imizes Fi under the constraint {xi < s,X2> X2) is {xi — s,X2 — ^2)- In this 
configuration Fi — —{k + ki2)s + ki2X2 which is strictly larger than zero since 
X2 > X2.0 

Lemma 3 F2{t) < and V2{t) < V20 <0. 

Proof: F2 = in the configuration {xi,X2) = {s,X2eq)- Since X2(t) > X2 > 
X2eq and Xi < s, the equations of motion guarantee F2{t) < 0. Then F2{t) < 
and V20 < imply V2{t) < V20 < O.n 

Lemma 4 xi{t) — s ^ vi{t) ^ 0. 
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Proof: For mass one to reach the stop with zero velocity it would have to be 
decelerating as it approached the stop. This never happens since by Lemma 2 
the mass always accelerates towards the stop.D 



Lemma 4 is important since it rules out the possibility of CD phases while 
X2{t) > X2 holds. In other words, every time mass one reaches the stop there will 
be a rebound. This fact will be assumed when deriving some of the lemmas that 
follow. 

In the lemmas that follow we initialize Vouto with viq. 



Lemma 5 



< \Vout„ 



Proof: Lemma 2 implies that mass one decelerates as it moves monotonically 
from s to some turning point a and then accelerates as it moves monotonically 
from a back to s, where the next rebound occurs. 

Now, V2{t) < f 20 < (Lemma 3) implies that, for each x e (a, s), Fi at x is 
smaller when mass one is going back to the stop than when it was coming out of 
the stop. Hence, the magnitude of the work done between s and a to bring the 
mass from li^oMtnl to rest is larger than the magnitude of the work done between 
a and s accelerating the mass from to Vin„_f_i ■ Therefore the kinetic energy of 
mass one just before the n+lth rebound is less than its kinetic energy just after 
the nth rebound. □ 



Lemma 6 



'^OUtn 



+1 



Proof: Combine lemma 5 and the rebound law Vout„ — —^'<^in„-^ 

2 

Lemma 7 xi{t) > xi where xi :— s — ^pl^in ■ 

Proof: By lemma 6 the maximum speed of mass one out of a rebound is 
boutol '■— 1^10 1- Now, assume mass one is acted throughout by the minimum 
possible force, F™*", as given in Lemma 2. This assumption yields the stated xi 
lower bound on xi{t).U 

Lemma 8 F^^"" < F2{t) where F^^"" := - (A; + k^) X20 + k^Xi. 

Proof: By lemmas 3 and 7, Xi{t) > Xi and X2{t) < X2o. The most nega- 
tive value of F2{t) occurs when m2 is at its rightmost position and mass one is 
at its leftmost position. Therefore, F^'", the value of F2 in the configuration 
{xi,X2) = {xi,X2o), is a valid lowerbound on F2{t). □ 

This ends the set of lemmas that hold while X2{t) > X2. 
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Lemma 9 . 

// 

1. The nth rebound occurred. 

2. The inequality Xi{t) > X2 holds throughout up to a time ^ after the nth 



1 

rebound 



Then 

The n+lst rebound is guaranteed to occur within that time interval, i.e., 

'-^^n — jpmin ■ 
^1 

Proof: This proof is by contradiction. Assume mass one has not yet reached 
the stop at the end of the time interval of length ^^^^/"L Then, by Lemma 2, 

< < Fi{t) throughout the time interval of length ^^^/"^ following the nth 

rebound. Therefore, at the end of this time interval, mass one must be moving 
again towards the stop with speed |fi| > I'^outnl- 

Now, let < s be the position of mass one at the end of this time interval. 
Noting that is strictly less than s, an argument parallel to the one used in 
Lemma 5 yields that |fi| at xl must be strictly less than l^outnl- Hence we get a 
contradiction. □ 

Lemma 10 . 

// 

X2(t) > X2 holds for t e [0, Er=o {^p^ 
Then 

the nr\-lth rebound occurs at a time t < J27=o ^ ^^m/n • 
Proof: Follows from Lemma 9 and Lemma 6 using induction. □ 
We are now ready to prove Theorem 1. 

Proof of Theorem 1: Let T be defined as the time it takes mass two to 
reach X2- We would like to bound T below by Ti. 

< Ti < T (3) 

By Lemma 8, this can be done by assuming F2{t) — F^*" while X2{t) > X2. Then 
Ti is just the (positive) solution to the kinematic equation 

X2 - X20 = V20T, + ^Fr^'T^ (4) 
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Note that Ti > as asserted because both F^^" and V20 are finite and X2 < X20. 

Now, suppose, instead of Vi(0) = Vio, we had picked an initiai condition 
v^{0) = vl';^''", where Vio < v^;^"" < 0. 

Looking at Lemma 7, one sees that our oid xi iower bound on xi{t) is stiff 
vafid. Likewise for the F^*" fower bound on F2{t) in Lemma 8. Therefore, Ti 
cafcufated in (4) is afso stiff a vahd fower bound on the time T"^"' it now takes 
mass two to reach X2, i.e., 

< Ti < T"""' (5) 



^ ' j = „ jz^ , it is cfear that we can pick a v 



new 
10 



vafue such that J2ilo y j < Ti is afso satisfied. Then, by Lemma 10, it 

foffows that mass one undergoes an infinite number of rebounds in the finite time 
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We have shown that for any X20 > £2, V20 < 0, we can find a f^g"' < 
such that the initiai condition (xi(0), Vi(0), 0:2(0), t;2(0)) = (s, v^q'", X20, 1^20) ieads 
to infinite rebounds in finite time. Now, if in any such {SjV^q"^ , X2o,V2o) initiaf 
condition we repface v'^q'^ by v^q, where v'^q'^ < v^q < 0, we stiff have an initiaf 
condition that feads to infinite rebounds in finite time. This is so because: 

• Ti is stiff a vafid fower bound for the time it now takes mass two to reach 

X2. 

/ O^i L,* I \ 

^00 



Fo J ^ ^i=0 y Fo 

Therefore we have a 3-D region in the 4-D phase space that feads to infinite 
rebounds in finite time. Noting that 

• This region is contained in the hypersurface {xi = s}. 

• Vi < for aff the points in this region. 

we concfude that orbits through this region evofve transversefy to it. Hence the 
region of phase space that feads to infinite rebounds in a finite time has nonzero 
measure in the phase space. □ 

2.3 Continuation of orbits that lead to infinite switching 

In order for an orbit that icads to infinite switching to be weff defined for aff time, 
we postuiate that orbit to continue from the point approached by the infinite 
rebounds. We caff this a reset of the orbit. 
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Definition 11 (Orbit Reset) Le^x(t;xo) be an orbit that leads to infinite switch- 
ing between COD dynamics and R dynamics as t — > t*. Assume thatx.{t; Xq) — 
X* as t — > t* . 

For t >t* we define x(t; xq) = x(t — t*;x*). We call x* a reset point. 

For this to be a valid procedure, one must show that a hmit point is indeed 
always approached under this form of infinite switching. 

Proposition 12 Let x(t,xo) be an orbit that leads to infinite switching between 
COD dynamics and R dynamics as t — > t*. Then, x(t,xo) — > x* as t — > t* , 
where x* = (s, 0, x*2, v^). 

Observe that, from the boundcdncss of energy, it follows that Vi{t) , X2{t) 

and V2{t) are bounded for t G [0, T). From this, in turn, it follows that F2{t) and 
Flit), except at rebounds, are also bounded for t e [0,r). 

Now we prove Proposition by proving a series of lemmas. 

Lemma 13 Let x(t,xo) be an orbit that leads to infinite rebounds as t — > t*. 
Then xi{t) — > s as t — > t*. 

Proof: Since there arc infinite rebounds ast^t*, one can find a time arbitrarily 
close to t* at which xi{t) = s. Together with the boundedness of this implies 
that Xi{t) — > s as t — > t*. 

Lemma 14 Let x(t, xo) be an orbit that leads to infinite rebounds as t — > t*. 
Then Vi{t) — > as t — >t*. 

Proof: First we prove that Lirrit^t* Inf > 0. Assume not. That is, assume 
that 3e > s.t. Vii,3i2 with ti < ^2 < t* and ^1(^2) < — e. Now, because -Fi(t) 
is bounded under COD, it is possible to bound away from the time it takes 
for mass one to go from a velocity equal or more negative than — e to rest. This 
bounds acts also as lower a bound on the time it takes for mass one to reach the 
stop again. But we can find a t^ arbitrarily close to t* satisfying Vi{t2) < — e. 
This contradicts the existence of infinite rebounds as t — t*. 

Now we prove that Limt^f Sup vi{t) < 0. Assuming otherwise, together with 
Xi{t) s as t ^ t* and the boundedness of Fi(t) under COD dynamics, implies 
that there exists an e > such that, arbitrarily close to t*, one can find a rebound 
with a velocity larger than or equal to e. Considering the velocities after those 
rebounds, one see that this violates Lirut^t* Inf vi{t) > O.D 

Lemma 15 Let x(t,xo) be an orbit that leads to infinite rebounds as t — > t* . 
Then there exist values X2 and such that X2{t) — > X2 as t — > t* and V2{t) — > 
V2 as t — > t* . 
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Proof: V2{t) has a limit since F2{t), i.e., the derivative of V2{t), is bounded for 
t e [0,t*). X2{t) has a hmit since V2{t), i.e., the derivative of X2{t), is bounded 

for t e [0,r).a 

The proof of Proposition 12 follows at once from lemmas 13, 14 and 15. 

We have therefore solved the problem of infinite switching between COD dy- 
namics and R dynamics in a finite time. Next, we would like to rule out any form 
of infinite switching in finite time involving CD dynamics segments. This will 
follow as a corollary to the proposition below. 

Proposition 16 It is possible to bound away from zero the time elapsed between 
any two CD dynamics segments. 

Proof: First note that F2{t) < for X2{t) E {x2eq, X2] and that this is vahd under 
all governing dynamics and regardless of the position of mass one. 

Now, from chart 2.1, wc know that at the end of a CD segment X2{t) = X2 
and V2{t) < 0. Therefore, after a CD segment, mass two will move monotonically 
from X2 to X2eg- Since V2{t) is bounded, we can bound away from zero, say by T, 
the time it takes mass two to go from X2 to X2eq- Note that this analysis still holds 
even if orbit resets due to infinite switching between COD and R dynamics occurs 
in the process. This is the case because an orbit reset maintains the continuity 
of X2{t) and V2{t). 

Now, the next CD segment can be joined either following a COD segment or 
following an orbit reset due to infinite switching between COD dynamics and R 
dynamics. In the first case the CD segment can be joined only at a point where 
mass two is at X2 again. In the case of an orbit reset at t* it can only be joined 
if limt^t*X2{t) > X2. In either case T is clearly a valid lower bound on the time 
elapsed between two CD segments. □ 

Corollary 17 Infinite switching in a finite time involving CD dynamics is not 
possible. 

In other words, infinite switching in finite time occurs only between COD 
dynamics and R dynamics. Definition 11 {Orbit Reset) is therefore all that is 
needed as far as defining continuation of orbits that undergo infinite switching 
between different dynamics in a finite time. 

The last step towards establishing forward existence and uniqueness is to show 
that there can be no infinite orbit resets themselves in a finite time. We start 
with a lemma about what are possible reset points. 

Lemma 18 Let x(t;xo) be an orbit that leads to infinite rebounds as t — > t*. 
Let x(i) — > x* as t — > t*. Then x* must be of the form {s, 0, X2, V2), where X2 is 
equal or greater than X2 . 
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Proof: That xl = s and t'jf = is shown in Proposition 12. Now, assume < X2. 
Then, X2{t) ^ X2 a.s t — > t* implies that it is possible to find a constant ci and 
a time ti before f* such that, for t G X2{t) < Ci < X2. In turn this means 

one can find a nonempty interval (c2,s) such that, for all t £ {ti,t*), when Xi{t) 
is in (c2, s), Fi is negative. 

Now consider any rebound occuring at a time t e {ti,t*). By the above 
analysis, after the rebound, mass one will have to leave the region (c2,s) before 
it can turn around and return for another rebound. Since Vi{t) is bounded, one 
can bound away from zero the time it takes for mass one to go from s to C2. This 
contradicts the existence of infinite rebounds as t — > t*. □ 

Corollary 19 After an orbit reset one of the following two types of segments is 
joined: 

1. A CD segment (this is the case if x% > xi or {x\ — xi and V2 > 0) ). 

2. A COD segment starting at a point of the form x* = (5,0,5:2,^2) where 
vl < 0. 

Proposition 20 It is not possible to have infinite orbit resets in a finite time. 

Proof: Each orbit reset is followed by a segment of one of the two types out- 
lined in the corollary above. Let us look at each possibility and its implications 
separately. Consider an orbit reset that leads to a COD segment starting at a 
point of the form x* = (5,0,5:2,^2), where V2 < 0. First note that F2(t) < for 
X2{t) e {x2eq,S;2\ and that this is valid under all governing dynamics and regard- 
less of the position of mass one. Together with X2 — X2 and V2 < 0, this allows us 
to conclude that mass two will be moving monotonically from X2 to X2eq- Since 
V2{t) is bounded, we can bound away from zero, say by T, the time it takes mass 
two to go from 5:2 to X2eq- As explained in the proof of proposition 16, this holds 
even if orbit resets occur in the meantime. Now, the next orbit reset leading to a 
segment of this type can only occur when the limit of 0:2 (t) is again 5:2. Therefore 
T is also a bound on the time elapsed before the next orbit reset which leads to 
a segment of this type. Thus, in a finite time, there can only be finitely many 
orbit resets that lead to segments of this type. 

Now, consider an orbit reset that leads to a CD segment. By proposition 16, 
the time elapsed between two CD segments is bounded away from zero. Therefore, 
so is the time elapsed between two orbit resets that lead to CD segments. Hence, 
in a finite time, there can only be finitely many orbit resets that lead to a CD 
segment. 

Putting the two cases together, one concludes that in a finite time there can 
only be finitely many orbit resets. □ 
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We conclude that the addition of postulate 11 to the mathematical model 
ensures forward existence and uniqueness of solutions for all time. 

Note that, however, backwards uniqueness is broken by orbit resets, since 
when an orbit undergoes a reset it joins a regular orbit of the system. In fact, 
orbit resets are the only way in which backwards uniqueness is broken in the 
system. 



In this second part of the paper we analyze the dynamics of the system. Rather 
than providing a long and exhaustive characterization, we highlight some of the 
more interesting features of the system. However, for completeness, we first state 
a general theorem concerning the forward limit set of the system. 

Theorem 21 The forward limit set of the system is the set of orbits that are gov- 
erned exclusively by coupled oscillator dynamics and contact dynamics, in other 
words, the set of orbits that do not go through any rebound. 

Since the theorem is quite plausible and the proof is long, wc ommit the proof. 
It relies on a LaSalle standard type of argument^. We just point out that, to use 
such an argument, invariance of the forward limit set must first be proved. Note 
that this is not immediate, because for a discontinuous system the forward limit 
set is not necessarily invariant, as the example below illustrates. 

Example 22 A discontinuous dynamical system with a non-invariant forward 
limit set. 



Clearly, the only point in the forward limit set of this system is the point 
X = 0. However, the forward limit set is not invariant since x = is not a fixed 
point of the dynamics. □ 



3.1 Periodic hybrid orbits 

In this section we study orbits which contain both CD segments and COD seg- 
ments but no rebounds, that is, orbits where mass one alternates finite periods 
against the stop with periods in which it oscillates together with mass two under 
coupled oscillator dynamics. Since mass one must avoid any rebounds, such orbits 
might seem unlikely. However they do exist. It turns out that all such orbits are 

^See Khalil [8]. 



3 The Dynamics of the System 
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periodic and consist of exactly two and only two segments, one governed by the 
contact dynamics, the other by the coupled oscillator dynamics. We call them 
2-phase periodic hybrid orbits. We start by showing why these 2-phase periodic 
hybrid orbits are the only possible orbits made of just CD and COD segments. 

Let us track such an hypothetical hybrid orbit, starting with a CD segment. 
First note that, following a CD segment, a COD segment can only be joined at 
a point in the set {s,0,X2,V2 < 0}. This is so, since these are the exiting points 
of the CD region (see chart 2.1). Let the COD segment be joined at the point 
(s, 0, 5;2, '^'2) ) where V2 is a negative number. Now, this point, when evolved 
under COD in backwards time, comes from the {xi > s} region. This means 
that, whether the COD orbit is periodic or quasiperiodic, if the system stays long 
enough under coupled oscillator dynamics, the orbit will return to the {xi > s} 
region. But we know this can never happen, so eventually a rebound or a CD 
segment must follow that COD segment. Assume a CD segment rather than a 
rebound follows the COD segment. This next CD segment can be joined only at 
a point on the entrance surface of the CD region/exit surface of the COD region. 
Prom chart 2.1, this must be a point of the form (s, 0, X2, ^2*), where V2* is a pos- 
itive number. Furthermore, because energy was conserved during the preceeding 
COD segment, we must have that V2* = —v^. Next, during the CD segment mass 
two oscillates as a harmonic oscillator for a fraction of a complete cycle. If the 
CD segment started at (s,0, £2, —'^2)) it must end at (s,0, X2,i'2) (see Figure 2). 
At this instant the orbit as exactly the same initial conditions as at the beginning 
of the description. Therefore everything simply repeats. We conclude that 

The periodic orbits described above are the only possible orbits that contain 
both Coupled Oscillator Dynamics and Contact Dynamics and no Rebounds at 
all. Conservation of energy precludes the possibility of more complicated orbits, 
periodic or otherwise. 

We now turn to the task of explicitly constructing 2-phase periodic hybrid 
solutions analytically. 

Let the solution to the COD segment of the hybrid orbit start at time i = 
and be expressed as 

xiit) = 3f?( Jie^-^i* + J2e^'^2i) 

X2{t) = 3?( Jie^'^i* - 126'''^^) ^ ' 

where ui — y/k, 102 — \/k-\- 2ki2 and Ji and J_2 are complex constants related to 
initial conditions by 
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X2(0) = K(Ji)-^U2) 



This is a segment of an hybrid orbit if and only if the following set of five 
equations and one inequality are satisfied for some integer n & N. 

s = ^Ui) + UU2) (8) 

- -UJ,QUl) - ^2^>U2) (9) 

X2 = 3?Ui) - ^{12) (10) 

V2{0) = -Ui^U,) + i02^{J2) (11) 



27r - 2 arctan ^ ^ 

27r - 2 arctan |J=4 + 27m ^2 



(12) 



27r — 2 arctan ^7j4 

x-i{t) < s is satisfied Vt e [O,icod], where tcod — (13) 

tcod being the duration of the COD phase of the hybrid orbit. 

(Note: the arctan range in these and all the subsequent formulas should be 
taken as [0, tt].) 

In the above set of equations, (8) ensures that mass one is initially at the 
stop, (9) ensures that its initial velocity is zero and (10) guarantees that mass 
two is initially at X2- Equation (11) simply relates ^2(0), the initial velocity of 
mass two, to the initial condition variables J_i and J_2. No particular requirement 
is imposed on the value of ^2(0). Equation (12) guarantees that xi{tcod) = xi{0) = 
s, vi{tcod) = vi{0) = 0, X2{tcod) = X2{0) = X2 and V2{tcod) = -V2{0). Finally, 
inequality (13) ensures that mass one stays to the left of the stop throughout the 
COD segment. The first five equations can be reduced to one by eliminating the 
complex variables J_i and J_2. The resulting equation is transcendental involving 
the constants k, ki2 and s and the required initial condition of mass two, ^2(0). 

27r— 2 arctan 



27r+2n7r-2 arctan -'='2 "2^0) Vfe+2fci2 , . 



27r-2 arctan- "'2^°> 



^ xi{t) < s : yt e [O,tcod], where tcod = 
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So, for any integer n for which one can find a value of 1^2 (0) that satisfies 
(14), one has a vahd COD segment of a 2-phase periodic hybrid solution. The 
constants J_i and J_2 in the analytical solution (6) of the COD segment are then 



■2 ~ 2k^ '■2x/fc+2fci2 



and this solution is valid for t e [O,tcod], the duration of the COD segment. 
The analytical solution to the CD segment of the orbit is 



f Xi(t) = s 

\ X2{t) = X2eq + ^{Ce''^^h--^*-*-''d)) 

where Ushm = ^/kTh2. Q = (^)s - X2eg + 
and this solution is valid for 

"2(0) 



(16) 



t e [tcod, tcod + ted], where ted 



27r - 2 arctan 



T. S 



^shm 



The constant C was chosen so that the boundary conditions 

1 X2{tcod) = X2 

\ V2{tcod) = -V2{0) 

between the COD and CD segments are met. 



(17) 



An example of a 2-phase periodic hybrid orbit in a system with spring con- 
stants k — 13, ki2 = 25 and the stop placed at s = 1.8 is shown in Figure 3. The 
dotted line represents the trajectory of mass one. The continuous line represents 

the trajectory of mass two. In this case ^2(0) = —17.47 satisfied the transcenden- 
tal equation in (14) for n=l. The times tcoci=l-289 and ted — 0.26, so the period 
of this hybrid orbit is tcod + ^cd = 1-549. 



3.2 Dynamics of action angle coordinates and simulations 

In this section we analyze the dynamics in the action-angle coordinates associ- 
ated with the COD. This allows a simple geometrical way of interpreting and 
visualizing the system dynamics. The analysis in this section concerns only the 
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Figure 3: A 2-phase periodic hybrid orbit. 



case in which the COD orbits are quasiperiodic. Assume first that there is no 
stop. Parametrize the tori in which COD orbits hve by (Ji, I2), the value the 
COD action variables take on those tori. A COD orbit expressed as in (6) lives 
in the torus where the action variables take the values Ii = \J_i\ and I2 = iJjzl- 

Definition 23 . 

Referencing orbits by the torus in which they live captures most of the in- 
formation, since quasiperiodic orbits in the same torus are qualitatively similar. 
Therefore we will now study the evolution of orbits solely in terms of of the evo- 
lution of their COD action coordinates. We call the space defined by the Ii and 
I2 action coordinates the action space. Naturally, without the stop, orbits never 
leave the torus they start in, i.e., they remain at a fixed point in the action space. 
Now we describe how the introduction of the stop changes this. The effect of 
introducing the stop is to divide the tori into three distinct groups: 

1. Tori satisfying Ii + I2 < s. 

These are tori that are a finite distance away from the {xi = s} surface. We 
call them interior tori. The orbits living in the interior tori are the ones 
where mass one stays always a finite nonzero distance away from the stop. 

2. Tori satisfying Ii + I2 = s. 

These are tori that are tangent to the {xi = s} surface. We call them grazing 
tori. Grazing tori are the tori that touch the hypersurface {xi = s} exactly 
at one point, of the form (s, 0, X2, 0), where X2 is between —s and s. As the 
name indicates, on orbits living in grazing tori mass one comes arbitrarily 
close to the stop, on some of those orbits actually touching it with zero 
speed. 
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3. Tori satisfying Ii + I2 > s. 

These are tori that are cut by the {xi = s} surface. 

This division is shown in Figure 4. The tori on the diagonal hne going from 
(0, s) to (s, 0) are the grazing tori. The tori to the right of the hne are cut by the 
stop, while the ones to the left of the line are the interior tori. 




Figure 4: Action space. 



An orbit that starts in a torus that is cut by the stop winds around that torus 
until it intersects the surface where the torus is cut. At that point, physically, a 
rebound occurs. In the action space, the orbit simply jumps to another torus at 
a lower energy level. It then starts winding around the new torus until it again 
reaches the surface where the torus is cut, at which point the orbit jumps again 
to a torus of lower energy. The jumps are never to one of the interior tori since 
interior tori form an invariant set. Also, because of quasiperiodicity, an orbit in 
a torus which is cut by the stop cannot remain under COD indefinitely without 
intersecting the {xi = s} surface. Orbits in these tori are therefore guaranteed 
to eventually jump to a lower energy torus. The exception to this is the 2-phase 
periodic hybrid orbit. The COD segments of 2-phase periodic hybrid orbits are 
in tori that are cut by the stop. However, when a 2-phase periodic hybrid orbit 
reaches the {xi = s} surface, it joins a CD segment instead of undergoing a re- 
bound. Then, when the CD segment is over, the exact same COD segment that 
lead to it is rejoined. The 2-phase periodic hybrid orbits are therefore these very 
special orbits, the only that manage to survive in tori that are cut by the stop. 

Figure 5 is a numerical simulation of the evolution of an orbit in the action 
space. The tori which the orbit visits have been connected by straight fines in 
order to allow one to discern the sequence in which the tori are visited. In this 
simulation the stop was placed at s = 4, the coefficient of restitution was c=0.8 
and the spring constants were A; = 5.3657 and A;i2 = 1.5682. The initial condition 
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of the orbit shown is (a;i(0), t'i(O), 0:2(0), ^2(0)) = (—6,0,6,0). As the orbit jumps 
from torus to torus, it asymptotically approaches the line of grazing tori. As it 
does so, the jumps become smaller and smaller. This can be seen in Figure 6, 
a zoom of Figure 5 near the grazing tori line. Predicting which grazing torus, if 
any, the orbit will specifically approach, remains an open research question. 




Figure 5: A trajectory in action space. 




Figure 6: Zoom of Figure 5 near the grazing tori line. 

Figure 7 shows a few more simulations that illustrate a range of the behaviors 
that can be observed. On each simulation exactly one parameter was changed 
from the parameters used in the simulation in figure 5. In the top left picture, a 
different initial condition, (a:i(0), fi(0), 2:2(0), ^2(0)) = (6,0,6,0), was used. On 
the top right picture the coefficient of friction was reduced to £ = 0.4. Hence, 
more energy is lost at each impact and within a small number of impacts the 
system comes very close to a grazing torus. On the bottom left picture, the stop 
was placed at s = 2. Interestingly, this orbit seems at first to be approaching a 
torus near (/i,/2) = (2,0) but then suddenly, with two or three larger jumps, 
moves to the region near {hjh) = (1.5,0.5). On the bottom right picture, the 
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spring constants were changed to A; = 2.5465 and ku — 3.9780. The trajectory 
in this picture clearly undergoes three qualitatively distinct phases. First there 

are jumps between far apart tori. Then, abruptly, the trajetory changes to very 
small jumps which gradually lead the system to the region next to (/i, I2) = (4, 0). 
Finally, once close to (/i, I2) = (4, 0), the trajectory slowly creeps up, parallel to 
the grazing tori line. 




li '1 

Figure 7: Trajectories under different parameters and initial conditions. 

We point out that this section's analysis is not pertinent when the COD orbits 
are periodic rather than quasiperiodic. In that case, there exist COD orbits that 
never intersect the {xi = s} surface even though they live in tori that are cut by 
the stop. 
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4 Conclusions 



In this paper we studied a coupled oscillator system with inelastic impacts. We 
started by showing that applying the standard model in which the system switches 
between the different types of governing dynamics results in a fundamental prob- 
lem, namely the existence of a nonzero measure set in the phase space where 
infinite rebounds occur in a finite time. This problem was solved by continuing 
such orbits from the limit point approached by the infinite rebounds. Then, look- 
ing at the system dynamics, we found a particularly interesting type of periodic 
behavior, the 2-phasc periodic hybrid orbits, which we described analytically. Fi- 
nally, we presented some numerical simulation results and a way of visualizing 
and interpreting them, based on action-angle coordinates. 

This paper highlighted two aspects of hybrid systems. First, it stressed the 
care that must be taken to ensure a well defined mathematical model for the 
hybrid system. Secondly, it demonstrated that hybrid systems, even elementary 
ones, can possess quite interesting and unusual dynamics. We hope this will serve 
as a useful step towards further theoretical studies of multi-degree of freedom 
hybrid systems. 
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